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ABSTRACT 

We use Chandra observations to estimate the accretion rate of hot gas onto the central supermassive black hole in four 

giant (of stellar mass M* ~ 10^^ — IO^^Mq) early-type galaxies located in the Virgo cluster. They are characterized by 

an extremely low radio luminosity, in the range L < 3 x 10^* — lO^'^ ergs~^ u,-i 

in these objects occurs at an extremely low rate, 0.2 — 3.7 x lO~^A/0 yr~^ 

accretion - jet power found for more powerful radio-galaxies. This confirms the dominant role of hot gas and of the 

galactic coronae in powering radio-loud active galactic nuclei across ~ 4 orders of magnitude in luminosity. A suggestive 

trend between jet power and location within the cluster also emerges. 

Key words, galaxies: active, galaxies: jets, galaxies: elliptical and lenticular, cD, galaxies: ISM 



^ Hz ^. We find that, accordingly, accretion 
and that they smoothly extend the relation 



1. Introduction 

The origin of the energetic processes occurring in active 
galactic nuclei (AGN) has been one of the most relevant 
astrophysical problems in the last decades. It is well es- 
tablished that the main driving mechanism is the release 
of gravitational energy of flows accreting on supermassive 
black holes (SMBH) for both radio-loud and radio-quiet 
AGN. While in the latter class the energy is mainly re- 
leased into electromagnetic radiation, in radio-loud AGN 
a substantial fraction of the accretion power is converted 
into the kinetic power of relativistic jets. What is still a 
matter of debate is the origin of the accreting material that 
fuels the AGN, how this is channeled onto the SMBH, and 
whether and how it is eventually accelerated into a rela- 
tivistic collimated outflow. 

X-ray observations of radio-loud galaxies recently indi- 
cated that in these objects spherical accretion of hot gas can 
provide the nece ssary level o f accre tion to power the active 
nucleus. Indeed, lAllen et all (|2006D found a strong correla- 
tion between the jet and the accretion powers in a sample 
of nine objects. The jet kinetic energy was estimated from 
the mechanic work needed to inflate the cavities observed 
in their X-ray images. The same X-ray data also provided 
information on the physical parameters of the hot inter- 
stellar medium (ISM), yielding estimates of the accretion 
power exploiting the Bondi model (1952). 

iBalmaverde et all (2008, here after BBC08) extended 
Allen's work considering a larger number of radio-galaxies. 
Because X-ray cavities are visible only in a minority of ob- 
jects, they took advantage of t he relati on between jet and 
radio-core power calibrated bv lHeinzet al. (2007i i. In this 
approximation, the jet kinetic power Pj is estimated di- 
rectly from the radio core luminosity, a quantity readily 
measurable in radio-galaxies. For a sample of 23 objects 
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with X-ray data of sufficient quality (including the sources 
analyzed in Allen et al.) they showed that the correlation 
between the jet and accretion power holds over '^ 3 decades, 
down to a jet power of ^ 10^^ ergs^^. Another result of 
this correlation is moreover that on average ^ 1% of the 
available accretion power is converted into jet kinetic en- 
ergy, independently of the AGN luminos i ty. Si milar results 
were obtained also by iHardcastle et all (|2007D following a 
different approach. 

Taking advantage of the recent surveys of early-type 
galaxies in the Virgo cl uster per formed in the optical 
(Cote et al. 2004) radio (jCaDetti e t al. 2009), and X-rays 
(jGallo et al.ll2OO80 we are in the position of testing the va- 
lidity of the correlation accretion/jet powers by including 
objects of even lower luminosities than those discussed in 
the previous studies. We then isolated four giant galaxies 
(of stellar mass M* ~ 10^^ — IO^^Mq) characterized by 
an extremely low radio luminosity (two of these sources 
are actually undetected in radio observations), in the range 
L < 3 X 10^^ — 10^^ erg s~^Hz~^, more than 5 orders of 
magnitude lower than the brightest objects considered by 
Allen et al. (2006 ) and BBC08. 

Furthermore, we have now a total of 10 radio-galaxies 
located in the Virgo cluster. It will then be possible to test 
whether in this environment the connection between the jet 
and accretion rates also holds and whether and how galax- 
ies show different properties, depending on their positions 
inside the cluster. 

In the following Sect. 2 we outline the main properties 
of the selected galaxies, in Sect. 3 we summarize how we 
estimated, through the X-ray and radio data, the values 
of the jet and accretion powers. In Sect. 4 we discuss our 
results and their implications for the sample considered by 
BBC08 with the addition of the four Virgo galaxies exam- 
ined here. The summary of our results are reported in Sect. 
5. In the Appendix we provide the details of the analysis of 
the Chandra X-ray data. 
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Table 1. Properties of the sample of core galaxies in the VCC 



Source 


Alt. names 


z 


D 


sub cl. 


d 


14 


Log Lcorc 


Log Mbh 


VCC 731 


NGC 4365 


0.00415 


23.3 


B 


6.2 


246 


< 25.84 


8.59 


VCC 798 


M 85, NGC 4382 


0.00243 


17.9 


A 


1.8 


578 


< 25.60 


7.93 


VCC 881 


M 86, NGC 4406 


-0.00081 


17.9 


A 


0.7 


1584 


25.94 


8.52 


VCC 1535 


NGC 4526 


0.00149 


16.5'' 


B 


> 1.1 


549 


27.03 


8.59 


VCC 763 


M 84, NGC 4374 


0.00354 


18.4 


A 


1.2 


247 


28.87 


9-OOir4'o 


VCC 1226 


M 49, NGC 4472 


0.00333 


17.1 


B dom. 


- 


- 


27.10 


8.78 


VCC 1316 


M 87, NGC 4486 


0.00436 


17.2 


A dom. 


- 


- 


29.99 


9.53t»;f, 


VCC 1632 


M 89, NGC 4552 


0.00113 


16.2 


A 


1.0 


967 


28.24 


8.54 


VCC 1939 


NGC 4636 


0.00313 


15.6'' 


B 


2.7 


59 


27.02 


8.16 


VCC 1978 


M 60, NGC 4649 


0.00373 


17.3 


A 


0.9 


190 


27.82 


9.30 in^ 



Column description: 1) name; 2) alternative name; 3) spectroscopic redshift from the NED database; 4) distance in Mpc (±0.5 
Mpc) from Mei et al. (2007); 5) sub-cluster to which each galaxy belongs, with the dominant galaxies marked; 6) and 7) three 
dimensional distance (Mpc) and velocity (km s~^) along the line of sight with respect to the dominant components of the sub- 
clusters Virgo A or B; 8) radio core luminosity ( ergs~^ Hz~^) from Capetti et al. (2009); for NGC 4636 we adopted the value 
of BBC08, adjusted at the distance assumed in this paper; 9) mass of the central SMBH in units of Mp derived fr o m the stellar 
velocity dispersion reported in the HyperLeda database or from direct measurements as listed bv lMarconi fc HuntI (|2003l ) (when 
not reported, we adopt an intrinsic error of ± 0.23 dex, related to the dispersion of the A4"bh - c* relation). 

^ For this object there is no distance measurement thus we adopted the average Virgo distance to estimate luminosities; accordingly 
its three dimensional distance d from the cluster is a lower limit, equal to the offset from M 4 9 in the plane of the sky. 
^ The distance of NGC 4636, with an assoc i ated e rror of ±0.9 Mpc, is fr om lTonrv et al.] l|2001f ). scaled by a factor 1.06, the average 
ratio between the distances of iTonrv et al.l (|200lf ) and lMei et al.l (|2007f ) for our sample. 



2. The sample 

We considered the sample observed by ICote et all (|2004l ) 
for their Hubble Space Telescope survey of the Virgo clus- 
ter (ACS/VCS). They selected a sub-set of 100 early- 
type galaxies from the Virgo cluster catalog (VCC, 
iBinggeli et al.l Il987l ) that consists of 2096 galaxies within 
this K, 140 deg^ region. These data provide us with 
a detailed analysis of their optical brightness profiles 
(jFerrareseet al.l(2006D. 



Balmaverde fc Capettil ()2006[ ) found that the nuclear 



multiwavelength properties of early-type galaxies are re- 
lated to their optical brightness profile. In particular, 
radio-loud AGN are found only in 'core-galaxies' (hereafter 
CoreG), i.e. objects whose brightness profile is better fit by 
a core-Sersic model compared to a pure 'power-law' Sersic 
model (hereafter SG ) extending to the innermost regions 
(|Truiillo et al.ll200l . 

We here focus on the properties of the sub-sample 
formed by the nine CoreG present in the Virgo samplqj- 
The VCC/CoreG are the brightest early-type galaxies in 
the sample (with the only pure Sersic interloper represented 
by 9th ranked VCC 1903). This was expected, considering 
the well kn own link between brightness profiles and lumi- 
nosity (e.g. iLauer et al.ll2007l ). Because the VCC is com- 
plete down to Br ~ 12 and all CoreG have B-r < 10.78, 
the VCC/CoreG sub-sample is also complete. These galax- 
ies can be well com pared with the radio-galaxies studied by 
lAllen et al.l (|2006[ ) and BBC08 considering the similarity of 
th eir host proper t ies. 

ICapetti et al.l ()2009D performed a radio survey of the 
Virgo Cluster, observing with the Very Large Array (VLA) 
in A configuration and at 8.4 GHz the 63 early-type galaxies 
covered by the ACS/VCS, limiting to those brighter with 



^ The low luminosity VCC 1250 is also classified as a core- 
galaxy but this classification is considered as marginal and 
should be interpreted with caution because the galaxy morphol- 
ogy is severely affected by dust patches and a bright nuclear 
source. 



Bt < 14.4. We here reproduce their Fig. 1, where the nu- 
clear radio luminosity is plotted against the stellar mass of 
the galaxy (see Fig, [p. In agreement with pr evious studies 
(e.g. lBest et"aIll2005riM"auch fc SadleiJl2007D . in the Virgo 
cluster the detection of a radio source is also highly favored 
in more massive galaxies. The nine VCC/CoreG have simi- 
lar masses, all within the range 10^^ — 10^^ ^^0- Conversely, 
they show a large scatter in their nuclear radio luminosity, 
over more than 4 orders of magnitude. Indeed, while seven 
are detected in radio, we could only set an upper limit for 
the remaining two. They differ even more in terms of tota l 
radio power (measured at 1.4 GHz, see lCapetti et al.|[2009[ ). 
spanning at least over a factor of > 2 x 10^. This confirms 
the idea that the relation between radio and optical lu- 
minosity can only be described in terms of a probability 
distribution. 

Below we explore the possibility that the energetic out- 
put in these galaxies is correlated with the accretion rate 
of the hot. X-ray emitting gas, in analogy with the results 
obtained on more powerful sources. Of particular interest 
are the four VCC/CoreG|j with the lowest radio luminos- 
ity (namely NGC 4365, M 85, M 86, and NGC 4526), not 
analyzed in previous studies, that reach levels as low as 
i-coro < 3 X 10^5 erg s-^Hz"!, a further factor of > 10 
lower than the faintest objects considered in BBC08. 

In Table [1] we list the main data of the four sources 
of interest. As we are also looking for possible connections 
between the accretion properties and the cluster environ- 
ment, we consider the relationships between the galaxies of 
our sample and the two main cluster's sub-structures (Virgo 
A and Virgo B, with M 87 and M 49 as dom inant compo- 
nents, respectively, see iBinggeli et al.) 119871 ). Accordingly 
we reported also the relative distances d and velocities V^ 
(along the line of sight) of the various galaxies with respect 
to the dominant galaxy of each sub-cluster. 



^ For easier reference we will use the more common Messier 
or NGC numbers for these sources. 
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Fig. 1. Radio core luminosity vs the galaxy stellar masses 
for the Virgo sample analyzed in Capetti et al. (2009). The 
four galaxies of interest are labeled with their Messier or 
NGC numbers. The radio detected galaxies are marked with 
large circles, separating them on the basis of the optical sur- 
face brightness profiles: filled circles are CoreG, while empty 
circles are pure SG. The undetected objects are marked 
with small squares (for clarity, we omit the downward point- 
ing arrows), empty for SG, filled for CoreG. 



At this stage, we also include in the sample the CoreG 
VCC 1939 (i.e. NGC 4636). This source is not covered 
by the ACS/VCS because it is reported by iBinggeli et al.l 
(|1987f) to be a possible cluster member, owing to its rela- 
tive large angular distance from ~ 8° (2.3 Mpc) from M 49. 
However, its membership appears to be rather secure, con- 
sidering its distance of 15.6±0.9 Mpc (compared to D=17.1 
Mpc of M 49) and its recession velocity of only 69 kms^ 
relative to M 49. In Table [1] we provide the data for easy 
reference also for the six VCC/CoreG already discussed in 
previous studies. 



3. Jet and accretion powers 

In this section we estimate the jet's kinetic energies and the 
accretion powers of the four CoreG with the lowest radio 
luminosity. The main assumptions and the details of the 
methods adopted to obtain Pb and Pj are outlined and 
discussed in BBC08 (and references therein). 



3.1. Jet power 

As already mentioned in the introduction, the jet kinetic 
power Pj can be estimated from the work done by the jets 
to inflate the ISM cavities observed in the X-ray images. 
This method can be used only in objects with high X-ray 
data quality and, o bviously, when such bubbles are visible. 
iHeinz et al.l (|2007[ ) showed that the radio core luminosity, 
Lcorc, can also be used to estimate Pj. Using a sample of 



radio-galaxies for which Pj could be measured from the 
presence of cavities, they calibrated a relation in the form 

Pj oc Lcoro that enables one to evaluate Pj from a sim- 
ple measurement of core radio power Lcorc, with a typical 
uncertainty (based on the rms of the relation) of ^ 0.4 dex. 

The four galaxies under scrutiny have rather different 
radio properties with respect to those considered by in 
previous studies. Indeed, while the galaxies considered in 
lAUen et al.l (|2006D and BBC08 have well-developed double 
radio structures, only NGC 4526 has an extended morphol- 
ogy indicative of a jet origin. The radio emission of M 86 
is instead unresolved, while it is undetected in NGC 4365 
and M 85. The very presence of jets in these last three 
sources cannot be taken for granted. We might refer more 
precisely to our analysis of the connection between accre- 
tion and jet power as a study of the link between accretion 
and the AGN radio emission. Nonetheless, because icorc 
and Pjct differ only for the scaling relation reported above, 
we prefer to maintain the jet power as a reference quantity, 
to facilitate comparison with previous works. 

This assumes that the Pj - Lcore relation holds for the 
sources discussed here, with a cor e power at le a st 10 times 
fainter than those considered by iHeinz et al.l (|2007l ). i.e. 
that their jets properties are similar to those of brighter 
radio-galaxies. In particular, there is the possibility that 
low-power jets are substantially slower, even non rela- 
tivistic. Implicitly adopting the average de-boosting factor 
for radio cores of more powerful sources for our sources, 
(5^ = 4 X r~^(2 — /3)~^, might cause an overestimate of 
the jet kinet ic power of up to a fac tor of ~ 3.6 (adopting 
r = 5 from iGiovannini et al.l l200lf) . Although this is not 
a negligible effect, it does not significantly affect our main 
conclusions. 

Because the estimates of Pj based on cavities can be 
considered as primary calibrators, while those derived from 
Lcore are secondary estimators, we preserved these more 
direct measurements for the nine objects of Allen et al. for 
the analysis that will follow in Sect. JD 



3.2. Accretion power 

The accretion power Pb = tjMbc^ is deduced assuming a 
steady, spherically symmetric m odel for accretion onto a 
gravitational body (|Bondilll952[ ). For the sake of simplicity 
we have assumed an efficiency rj = 1] the accretion rate, 
Mb, is defined at the Bondi radius tb = 2GMbh/cI (c^ is 
the sound speed, Mbh the mass of the SMBH) as 



Mb oc CsPbTb oc ubMbhTb 



-1-3/2 



(1) 



where pB, riB, and Tb are the mass and number densities, 
and the temperature at vb, respectively. 

The mass of the SMBH has been estimated from the 
stellar velocity dispersion, availa ble for all sources frorn . the 
HyperLeda database, using the iTremaine et al.l ()2002f) re- 
lation. In three CoreG of the Virgo sample, the mass of 
the SMBH has be en measured directly an d we adopted the 
value reported bv lMarconifc HuntI (|2003D . 

The temperature and density were deduced from the 
properties of the hot ISM through Chandra X-ray obser- 
vations. The data analysis and the results of these obser- 
vations are presented in detail in the appendix, where we 
derive the deprojected profiles of temperature and density 
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Table 2. ISM properties, accretion rates, jet and accretion powers. 



Name 


TB 


a 


Tb 


Ub 


Mb 


logPs 


log P, 


NGC 4365 
M 85 
M 86 
NGC 4526 


29 
5.5 
17 
19 


1 0+U.3 
n 8+0-3 

0.8_o2 
0.4«-i 


0.44 ± 0.08 
0.51 ± 0.23 
0.65 ± 0.02 
0.68 ± 0.21 


7+1-1 

"•'-0.4 

7+1-° 

"•'-0.6 

o.2«:^ 


4.6 
0.2 
0.5 
1.1 


1.41 ±0.38 
0.05 ±0.61 
0.45 ± 0.54 
0.80 ±0.62 


< -1.68 

< -1.83 
-1.61 
-0.83 


M 84 
M49 
M 87 
M89 

NGC 4636 
M 60 


47 
28 
142 
17 
9 
83 


n c:c:+0-iU 
"•^0-0.18 

o.oo«:}« 

r, oo+O-OS 

1 n+°-i 

J^-U_0.2 


0.71 ± 0.05 
0.90 ± 0.16 
0.80 ± 0.01 
0.67 ± 0.09 
0.54 ± 0.11 
0.80 ± 0.02 


o.9t;!- 
o.8t°:^3 

17+0-04 

"•^'-0.03 
1 q+3.5 
^•'J-0.8 

4 +0-3 

U-* -0.2 
"•^^0-0.04 


6.6 
12.0 
24.0 
4.4 
0.34 
11.0 


1 t^s+o-^a 

^^•O°-0.30 
1 OO+0.23 
^•»J-0.25 
r. 1O+0-40 
^•-'-'^-0.28 
1 qo+0.21 
^-•■J°-0.22 

0.30«-^l 

i.79«;r4 


U.O^_Q 29 

-0.80l»;io 

1 ro + O.lO 

-0.27 



Column description: 1) name; 2) Bondi radius in pc; 3) slope of the density profile; 4) temperature at the Bondi radius in keV; 5) 
number density at the Bondi radius in cm~'^; 6) Bondi accretion rate in 0.001 x Mq yr~^; 7) accretion power in 10'*'^ erg s~ ^; 8) jet 



power in 10*^ erg s ^ (when not reported a fixed error of 0.4 dex has been assumed, the other values are from lAUen et al.ll200' 
Table 3. Results of the statistical and correlation analysis. 



X 


Y 


N 


r 


P 


Pp 


a 


b 


rms 


pV 
^B 


P| 


10 


0.94 


0.79 


0.018 


1.17 ± 0.16 


-2.13 ± 0.23 


0.25 


Pb 


Pj 


27 


0.81 


0.79 


10-* 


1.05 ± 0.13 


-1.96 ± 0.27 


0.51 


pV,* 

^b 


pr 


10 


0.93 


0.79 


0.03 


0.91 ± 0.17 


-1.38 ± 0.16 


0.29 


n 


pj 


27 


0.62 


0.64 


0.002 


1.27 ± 0.21 


-1.24 ± 0.27 


0.69 


M^ff 


pi 


10 


0.86 


0.80 


0.016 


1.71 ± 0.36 


-15.6 ± 3.1 


0.37 


Mbh 


pj 


27 


0.59 


0.64 


1.2 X 10"^ 


2.20 ± 0.41 


-19.5 ± 3.6 


0.73 


mkpc 


pV 
^B 


10 


- 


0.49 


0.14 


- 


- 


- 


nikpc 


Pb 


27 


0.42 


0.39 


0.05 


1.67 ± 0.38 


4.10 ± 0.52 


0.80 



Column description: 1) dependent variable, and 2) independent variable (superscripts V indicate that only the Virgo sample is 
considered); 3) number of objects; 4) linear correlation coefficient; 5) Spearman rank correlation; 6) probability of obtaining the 
observed (or a higher) value of p under the null hypothesis of no correlation between the two variables; 7) slope, and 8) intercept of 
the correlation (Y = a X ± b; the two upper limits on Pj for NGC 4365 and M 85 are not considered); 9) scatter of the correlation. 
* The density ninn measured at the innermost annulus (see text) is adopted to estimate Pb- 
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Fig. 2. Jet power Pj vs the accretion power Pb (left) and the mas s of the SMBH Mbh {right) for the sample of 27 objects, 
consisting of the 23 objects from BBC08 and lAUen et al.l ()2006l) with the addition of the four new measurements. The 
10 sources belonging to the Virgo cluster are marked with red circles, while full circles indicate the four objects analyzed 
here. The solid lines represent the best-fit regressions for the whole sample. 
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across the coronae for the four VCC sources. The tempera- 
ture profiles across the coronal regions are essentially con- 
stant, with typical values of T « 0.4 — 0.7 keV. For the 
calculation of the Bondi accretion rate we assumed that 
Tb is equal to the temperature at the innermost shell. The 
value of ub was extrapolated down to r^ assuming a rela- 
tion n{r) (X r~", with a deduced from a power-law fit to 
the density profile of the ISM (see the appendix for details) . 
The values of these quantities, with the derived accre- 
tion and jet powers, are reported in Tabled which also in- 
cludes the data for the other Virgo objects from BBCO80 
We remark that for our four VCC objects the quality of the 
data is such that we were able to derive an estimate of Mb 
with a typical uncertainty of a factor of 3. Besides these 
purely statistical uncertainties, systematic effects can also 
play a role. For example we will discuss in the next Section 
the effects of the assumption of the power-law dependence 
of the density on radius. 

4. Discussion 

4.1. Accretion versus jet power 

In the left panel of Fig. [5] we compare the jet and the 
accretion power, i.e. Pj and Pb- The four radio- faint 
VCC/CoreG are all located in the bottom left corner of this 
diagram. This indicates that their low radio luminosity is 
associated with a low level of accretion and with values con- 
sistent with the correlation between Pj and Pb derived for 
larger powers. The results of the statistical analysis are re- 
ported in Table 3: the parameters of the linear correlations 
are deduced from the bisectrix of the fits interchanging the 
two related quantities (excludi ng upper liinits) a s indepen- 
dent and dependent variables (jlsobe et al.lll990[ ). 

The inclusion of these four objects only marginally al- 
ters the statistical parameters describing the Pj vs. Pb re- 
lation with respect to the results of BBC08 (see our Table 
[3] in comparison with their Table 3), yielding 

log Pj,43 = (1.05 ± 0.13) logPB,43 - (1.96 ± 0.27), 

where both powers are given in units of 10'*'^ erg s^^. 

Focusing on the sub-sample of 10 Virgo galaxies, Pb 
and Pj are also well correlated, with Spearman rank corre- 
lation p = 0.79, corresponding to a probability P^ =0.018 
of obtaining the observed (or a larger) value of p under 
the null hypothesis of no correlation. Slope and intercept 
agree within the errors with those found for the full sam- 
ple. Because these sources represent the majority of the 
low power sources, this indicates that the results of BBC08 
also hold for very low accretion rates and jet powers, and 
in particular that the conversion efficiency from accretion 
to radio or jet power does not significantly vary over ~ 4 
orders of magnitude. 

We also tested how critical our assumption is that the 
density ns can be extrapolated using a relation in the form 
n ex r"" down to the Bondi radius. Indeed, even though 
the density profiles in our sources appear to monotonically 
increase toward t he center (with the exception of M 87, see 
lAUen et al.ll2006l) in the inner regions which are not probed 
by the Chandra data, th ey might flatten (as pointed out 
bv iHardcastle et al]|2007D . We then explored the effects of 
adopting as Ub the value of the density at the innermost 



1) 



o 




A few misprints in their Table 2 have been corrected. 



''°g ri| ^p^ (in cm 3) 

Fig. 3. Pb vs nikpc, the density of the ISM at 1 kpc. 
Symbols as in Fig. [21 i.e. red circles represent the 10 sources 
belonging to the Virgo cluster, with full circles indicating 
the four new measurements. 



annulus of the profile (nin„), which can be considered as a 
strict lower limit to the actual density at r^. Despite the 
very large range in the value of the innermost radius ri^n 
compared to the Bondi radius, renormalizing the value of 
Pb with Ub — T^inn, the correlation Pj — Pb holds (with 
P^ and Pp = 0.03 and 0.002, respectively; see Tabled. 
The quasi linear dependence is maintained, although with 
a slightly higher scatter, and with an intercept (and thus a 
conversion efficiency) increased by a factor of ~ 5. 

Following BBC08, we looked for possible relationships 
between the jet power and the single quantities defining 
Pb through Eq. (1). No correlation was found between Pj 
and the temperature at the Bondi radius. A correlation 
Pj — fT-B emerges only when the whole sample is considered 
{Pp = 0.03) but with a large scatter (rms = 0.81), an indi- 
cation that the highest jet powers are in general found in 
the densest (at their Bondi radius) coronae. A stronger re- 
lationship of Pj with Mbh is conversely present as shown 
by the right panel of Fig. [2] and in Table |3l The slope of 
the fit (a « 2.2 ± 0.4 for the fufl sample and a w 1.7 ± 0.4 
for the Virgo sub-sample) agrees with that predicted by the 
dependence of the accretion rate from Mgj^. We must re- 
mark however that the scatter of the correlation Pj - Mbh 
is larger than for the Pj - Pb relation. From these argu- 
ments it is clear that the mass of the SMBH is the dominant 
parameter in driving the accretion process and powering 
the jet, but that the density at the Bondi radius plays an 
important role, reducing the scatter of the correlation. 



4.2. Pb 3nd large scale properties of the ISM 

iBalmaverde et all (J200 8^ found that the different levels of 
nuclear activity are related to global differences in the struc- 
ture of the galactic hot coronae. Indeed, a relation links the 
jet power with the corona X-ray surface brightness, albeit 
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with a large scatter. This suggests that a substantial vari- 
ation in the jet power must be accompanied by a global 
change in its ISM properties. As the temperature is quite 
uniform across the ISM, the density in the external region 
is the most relevant parameter that could be related to the 
accretion process occurring in the core. We here re-examine 
this issue, but instead of considering the surface brightness 
of the corona, we measured its de-projected density at a 
fixed radius of 1 kpc, nikpc, ~ 1 - 2 orders of magnitude 
larger than the Bondi radius, interpolating the neighboring 
points in the density profile. In Fig. [3] we compare nikpc 
with Pb- 

For the complete sample of 27 galaxies, the correlation 
analysis between Pb and nikpc provides Pp — 0.05 and 
a ~ 1.7±0.4|jThe residuals have a rms of 0.8 dex. A similar 
relation is present also between Pj and nikpc- The inclusion 
of the four low-power VCC galaxies is crucial to unveil the 
connection between these quantities, because it extends the 
range covered by both variables, providing the necessary 
leverage to the data. We thus confirm that the level of AGN 
activity is connected with the overall properties of the hot 
corona. 

The rather loose relation between accretion rate and 
coronal density is not surprising. Indeed, assuming typical 
sound speeds of '^ 500 km s""'^, the outer (1 kpc) and inner 
regions {tb ~ 10—100 pc) of the ISM are causally connected 
on dynamical timescales on the order of a few 10^ years. The 
large scatter observed can be related to events occurring 
inside the ISM on shorter timescales. These could be of 
external origin, e.g. transient changes due to mergers or 
interactions with the ISM, but the accretion flow can also 
be affected 'from inside', e.g. due to a variable release of jet 
kinetic energy into the corona. 

4.3. Accretion and the structure of the Virgo Cluster 

We now look for possible relationships between the jet 
power (strictly connected to the nuclear activity) of the 
Virgo CoreG and their location within the cluster, and in 
particular with respect to the dominant galaxies of Virgo 
A and B (see Table [T|). In Fig. [4] we compare Pj with (left 
panel) the three dimensional distance of each galaxy with 
respect to its relative sub-cluster dominant component and 
(right panel) also in terms of difference of recession velocity. 

A suggestive trend between distance from the dominant 
galaxy of each subcluster and jet power emerges; galaxies 
located closer to the center of each sub-cluster have a larger 
jet power than those at larger distances. Due to the large 
uncertainties and to the low number of objects this does not 
anyway correspond to a statistically significant correlation 
between Pj and d. An analogous trend appears when con- 
sidering the absolute value of the relative radial velocities, 
in the sense that high jet powers are not found in objects 
with the higher velocities. Clearly this analysis is limited 
because we are only able to measure one of the velocity 
components. 

The clearer picture is derived by looking at distance and 
velocity simultaneously. The case of M 86 is particularly in- 
structive. It is located at ^ 1 Mpc from M 87, at a distance 
similar to that of M 84, M 89, and M 60. With respect to 
these sources, M 86 has a jet power '-^6 — 60 times lower 



* No convincing correlation is found though considering only 
the 10 Virgo galaxies, see Table 3. 



(and a similarly lower accretion rate). This could be due to 
different dynamical effects in these galaxies. Indeed, M 86 
has the highest radial relative velocity with respect to M 87 
(Vz « 1600 km s~^, being actually blue-shifted with respect 
to the galaxy) and its Chandra image shows a long tail in 
the NW direction (jRandall et al.ll2008l ). a clear evidence of 
ram pressure stripping. The coronal gas content is progres- 
sively depleted at each passage in the central regions of the 
cluster, possibly causing the lower level of activity. 

From the available data it is difficult to separate these 
environmental effects from a possible mass segregation that 
would cause the more massive galaxies to be preferably 
found closer to the cluster's center. These would have corre- 
spondingly higher SMBH masses that due to the important 
role played by Mbh might be associated with objects with 
higher accretion and higher jet power. 



5. Summary and conclusions 

We have extended previous studies on the connection be- 
tween accretion of hot gas and jet power in radio-loud AGN 
with a sample of four objects belonging to the Virgo clus- 
ter, characterized by extremely low radio luminosity. Two 
of them are actually undetected in 8.4 GHz VLA observa- 
tions, corresponding to a limit to their radio luminosity of 
icorc ^ 3 X 10^^ erg s~^. These four galaxies are partic- 
ularly interesting because they have similar properties to 
those typical of the hosts of radio-loud AGN from the point 
of view of their brightness profiles (they are giant early type 
galaxies with a flat nuclear stellar core) and are associated 
with black holes of large mass (Mbh ~ 10® - lO^-^M©). 
Nonetheless, these objects are > 30, 000 times fainter than 
the radio core of M 87 (the brightest VCC galaxy) and 
> 200, 000 times fainter in total radio power. Among the 
various alternatives that can be proposed we here explore 
the possibility that this is due to an extremely low accretion 
rate. 

We used Chandra observations of these sources to esti- 
mate the power available from accretion of hot gas, Mb, 
in a spherical approximation. This was obtained by de- 
projecting the temperature and density profiles of each 
source. Because the diffuse emission in these targets is on 
average much fainter than in those considered in previous 
studies, we considered and removed the contamination from 
stellar sources to isolate the genuine ISM emission. The 
quality of the data is such that we could estimate Mb with 
a typical uncertainty of a factor of 3. 

Furthermore, we now have a complete sample formed 
by all 10 core-galaxies located in the Virgo cluster. It has 
then been possible to verify whether the connection be- 
tween jet and accretion rates holds also in this environment 
and whether galaxies show different properties, depending 
on their positions inside the cluster. 

Our results can be summarized as follows. 

1. The estimated accretion rates of the four radio- faint 
VCC galaxies are among the lowest of the sample consid- 
ered. This indicates that their low radio luminosity is most 
likely due to a correspondingly low level of accretion. 

2. The correlation between the accretion and the jet powers 
spans over 4 decades down to very low values of Pb and Pj, 
with an approximately constant conversion rate of Pj ~ 
0.01 X Pb. This even holds relaxing our assumption that 
the density profile can be extrapolated with a power-law 
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from the Galaxy cannot be measured using the brightness fluctuation method (due to a prominent dust lane) we only 
have a lower limit for d. NGC 4365 is instead indicated as an arrow with its large distance (6.2 Mpc) reported explicitely 
for easier visualization. 



to estimate its value at the Bondi radius. The mass of the 
SMBH is the main driving parameter of the correlation. 

3. The correlation between Pg and Pj is also found for the 
subsample of the 10 Virgo sources, with parameters similar 
to those describing the full sample. 

4. A rather loose relation links the accretion rate and coro- 
nal density at 1 kpc, suggesting that a substantial variation 
in the jet power must be accompanied by a global change 
in its ISM properties. 

5. A suggestive trend (that however does not correspond 
to a statistically significant correlation) is present between 
the location (and relative velocity) of a source with respect 
to the cluster's center and its jet power. In one case, this 
is clearly due to ram pressure stripping of hot coronal gas 
from a galaxy moving at high speed with respect to the 
cluster. 

The main findings of this work are additional support 
for a linear relation between the power available from ac- 
cretion of hot gas and the jets kinetic power and that the 
efficiency of this process is apparently independent of the 
level of luminosity, of the black hole mass, and of the envi- 
ronment. Although we are still far from understanding the 
dynamics of the accreting gas close to the SMBH, these re- 
sults represent robust constraints on the processes of jets 
formation. 

With the current instrumentation it is unlikely that a 
much larger range of AGN luminosity can be probed, con- 
sidering (at the low end) the depth of the available radio 
observations and (at the high end) the limit imposed by the 
spatial resolution of the X-ray images in more distant and 
powerful AGN. 

It would instead be important to detect the two galax- 
ies still elusive in our VLA survey in deeper radio images 



(currently within reach thanks to the Extended VLA) to 
confirm (or disprove) that they are genuine active galaxies, 
even though of very low luminosity. They can be very use- 
ful to probe the properties of accretion and jet production 
at extreme low powers, in particular NGC 4365 which is 
already a marginal outlier in the Pj — Pb diagram. 

In the objects considered until now the accretion pro- 
cess is highly radiatively inefficient and the electromagnetic 
output from the nucleus is negligible with respect to what 
emerges as kinetic power (BBC08). The situation is rather 
different for radio-quiet AGN, which are brighter by a fac- 
tor ^ 100 than r adio-loud AGN in th e X-ray band at equal 
radio power (e.g. iPanessa et al.ll2007n . It would be of great 
interest to explore whether, in analogy to the Pj — Pb re- 
lation in radio-loud sources, a correlation between the ac- 
cretion power and radiative output (e.g. a link between Pb 
and Lx) exists in radio-quiet objects. An initial attempt 
in this direction was made by ISoria et al.l (2006a, b) ana- 
lyzing a sample of eight galaxies selected for their low X- 
ray luminosity. No correlation was found, but these objects 
span only a factor of 10 in Pb and Lx- An extension to 
radio-quiet sources covering a much larger range of X-ray 
luminosity is clearly needed to firmly test this point. 
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Appendix A: X-Ray data and ISM properties of the 
four radio-faint VCC galaxies 

The Chandra Archive information for the four VCC galax- 
ies analyzed here are reported in Table lA.ll For the X-ray 
data reduction and analysis we followed the same proce- 
dures and software packages as in BBC08. To account for 
the possible contamination from the active nucleus, we ex- 
cluded the innermost region from the analysis, adopting a 
minimum radius of I'.'S. This is rather conservative, because 
a faint X-ray central point source is seen only in NGC 4365 
(|Sivakoffet al.ll2nn8h . 

Since the diffuse emission in these targets is on aver- 
age much fainter than in those discussed in BBC08, the 
contamination from stellar sources must be taken into ac- 
count. These sources mainly belong to two classes: low-mass 
X-ray binaries (LMXB) and stars of old populations with 
hot coronae, mainly cataclysmic variables (CVs) and active 
binaries (ABs). 

The brightest LMXB appear in the X-ray images as 
discrete sources superimposed to the diffuse coronal emis- 
sion: they have been localized and removed from the event 
files through the wavedetect command. However, less lu- 
minous, unresolved LMXB can still have a non negligible 
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Fig. A.l. Spectral fit to the spectrum of NGC 4365 adopt- 
ing a thermal (Mekal) plus a continuum power-law com- 
posite model, representative of the ISM diffuse and LMXB 
emission, respectively. 

contribution to the extended luminosity. As their emission 
is well described by a power-law model, they can be sepa- 
rated from the hot gas thermal emission through a spectral 
decomposition. 

For this purpose we performed a two-components spec- 
tral fit, Mekal + Pow Law, adopting the galactic values for 
the column densities Nh- For the thermal model the met- 
alicity is fixed to half the solar abundances, while for the 
power-law the spectral index has been fixed to the value 
found from the discrete sources for each galaxy. The result- 
ing power-law photon indices are in the range F ga 1.3— 1.6, 
i n goo d agreement with the measurements of llrwin et al.l 
(|2003D . We repeated this procedure considering four annuli 
for each galaxy. In three objects the integrated contribution 
of the unresolved LMXB is at a similar level of the thermal 
emission from the ISM, while it is negligible in NGC 4526. 
As an example we show in Fig. lA. li the results obtained for 
NGC 4365. 

The second contaminating component originates from 
the stellar coronae of CVs and ABs; their contribution can- 
not be separated through a spectral analysis as for LMXB, 
because thei r thermal X-ray emissio n is similar to that of 
the ISM fe.g. lRevnivtsev et al.l2007p. We explo it conversely 
the result obtained bv iRevnivtsev et al.l (|2008[ ): they show 
that the old stars contribution to the X-ray emission of 
early-type galaxies (per unit stellar luminosity in the K 
band) displays a small scatter (less than a factor of 2) 
around the value they measured for NGC 3379, 
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We then retrieved the K band surface brightness pro- 
files of the four galaxies of interest from the 2MASS Large 
Galaxy Atlas. Assuming that the IR-to-X-ray scaling rela- 
tionship also holds in our sources, we compared them with 
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Fig. A. 3. Deprojected temperature (left panels) and density (right panels) profiles of the four VCC galaxies. The two 
dotted lines in each figure are the linear fits on the innermost data points, used for the extrapolation of the density to 
the Bondi radius. g 
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Fig. A. 2. X-ray brightness profile of NGC 4365, as ob- 
tained from the Chandra observations (upper curve). The 
lower curve is an estimate of the contribution of X-ray emis- 
sion from old active stars calculated scaling the K band 
surface brightness profile by using Eq. A.l. 



the X-ray brightness profiles obtained from our analysis. 
In Fig. IA.2l we show as an example the case of NGC 4365. 
This is the galaxy where we measured the largest (possible) 
contribution from old stars that amounts to ^ 30 % of the 
total counts. We conclude that the diffuse soft X-ray light 
is dominated by the genuine emission from the galaxies hot 
coronae. Nonetheless, we corrected the observed profiles for 
this marginal contamination. 

After removing the contribution from stellar sources, we 
followed the same procedure to derive the deprojected pro- 
files of temperature and density as in BBC08, to which we 
refer the reader for further details. Very briefly, to mea- 
sure the dependence of temperature with radius, we de- 
project the spectra assuming spherical symmetry using the 
PRO JECT model in XSPEC, considering four annuh. The 
results are graphically reported in Fig. lA.3l Concerning the 
density, we first of all determined the X-ray surface bright- 
ness profile (SBP) extracting the counts in a series of cir- 
cular annuli. The next step is to de-project the observed 
SBP, i.e. obtaining the number of counts emitted per unit 
volume as a function of radius. Assuming spherical symme- 
try, the count contribution provided by each spherical shell 
to the inner on e s were determined following the calculation 
by iKriss et al.l ()1983[ ). From the deprojected SBP we ob- 
tained n{r) assuming thermal emission and inverting the 
normalization coefficient of the spectral model. 

As the Bondi radius is well inside the innermost annu- 
lus (by a factor between 3 and 20), the value of the density 
was extrapolated down to rg through a fit across the pro- 
file, assuming a power-law in the form n(r) oc r~". To limit 
systematic errors related to the arbitrary choice of the range 
of radii to be included in the fit, we performed the analysis 
twice by using four and six of the innermost density points. 
For each case we estimated the uncertainty related to the 
accuracy of the parameters of the fit describing the den- 
sity behavior. We forced the density at rs to be higher (or 



equal) than its measured value at the innermost annulus. 
We finally adopted for ub the average of the values found 
in these two extrapolations and for its uncertainty the full 
range given by the overlap of the two individual error bars. 
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